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Abstract. We consider the Stackel transform, also known as the coupling-constant metamorphosis, 
which under certain conditions turns a Hamiltonian dynamical system into another such system 
and preserves the Liouville integrability. We show that the corresponding transformation for the 
equations of motion is nothing but the reciprocal transformation of a special form and we investigate 
the properties of this transformation. This result is further applied for the study of the fc-hole 
deformations of the Benenti systems or more general seed systems. 



Introduction 

The Stackel transform also known as the coupling-constant metamorphosis (cf. e.g. also 
|25j). is a powerful tool for producing new Liouville integrable systems from the known ones. This 
is essentially a transformation that maps an n-tuple of functions in involution on a 2ra-dimensional 
Poisson manifold into another n-tuple of functions on the same manifold, and these n new functions 
are again in involution. In the present paper we show that the corresponding transformations for 
equations of motion are nothing but reciprocal transformations. We also study the properties and 
present some applications of the latter. 

The significance of reciprocal transformations in the theory of integrable nonlinear partial dif- 
ferential equations is well recognized. These transformations were intensively used in the theory 
of dispersionless systems as well as the theory of soliton systems (see e.g. |2TH 12*2] and references 
therein). However, the role of the reciprocal transformations in the theory of finite-dimensional dy- 
namical systems is far from being fully explored, and the goal of the present paper is to contribute 
to such an exploration by developing the theory of reciprocal transformations for Liouville integrable 
Hamiltonian systems. To the best of our knowledge, such transformations first appeared in the paper 
[Li] by Hietarinta et al., where the concept of the coupling-constant metamorphosis, or the Stackel 
transform [5] (cf. also |2S1 for even more general transformations in the action-angle variables and 
[03 EH El EH] for more recent developments), was introduced. The reciprocal transformation ap- 
peared in this context as a transformation expressing the time (evolution parameter) for the target 
system through that of the source system [T3] , but the question of whether it sends the (solutions of) 
equations of motion for the source system into those of the target system was not addressed in [Pl| . 

In fact, as we show below, this transformation, when applied to the equations of motion of the 
source system, in general does not yield the equations of motion for the target system, unless we 
restrict the equations of motion onto the level surfaces of the corresponding Hamiltonians, see Propo- 
sitions El and El below for details. 

Even more broadly, we show that two Liouville integrable systems related by an appropriate Stackel 
transform for the constants of motion are related by the reciprocal transformation for the equations 
of motion restricted to appropriate Lagrangian submanifolds (see e.g. Ch.3 of ^H] and references 
therein for more details on the latter). 
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Moreover, we present a multitime extension of the original reciprocal transformation from Hietar- 
inta et al. [T3|, and study the applications of this extended transformation to the integration of 
equations of motion in the Hamilton-Jacobi formalism using the separation of variables (cf. jU]). 

As a byproduct, we present reciprocal transformations for a large class of dispersionless, weakly 
nonlinear hydrodynamic-type systems, the so-called Killing systems || that are intimately related 
to the Stackel-separable systems O EH E| . 

In the second part of the paper we consider the relations among classical Liouville integrable Stackel 
systems on 2n-dimensional phase space. In [7j infinitely many classes of the Stackel systems related 
to the so-called seed class, namely, the fc-hole deformations of the latter, were constructed. Here we 
show that any fc-hole deformation consists of a sequence of elementary deformations (one-hole defor- 
mations). These elementary deformations are nothing but particular cases of the Stackel transforms 
considered in the first section of the present paper. Hence the equations of motion for infinitely many 
classes of Stackel systems are related, upon restriction onto appropriate Lagrangian submanifolds, 
to the equations of motion for the systems from the seed class (which is a natural generalization of 
the Benenti class for the classical Stackel systems) by a sequence of reciprocal transformations. 

The significance of this result stems from the fact that the overwhelming majority of known 
today Liouville integrable natural dynamical systems that admit orthogonal separation of variables 
belong to the seed class. The /c-hole deformations of such systems again are Liouville integrable 
natural dynamical systems and admit orthogonal separation of variables, but the corresponding 
separation curves no longer are of the seed type. Our results make it possible to understand how the 
corresponding dynamics is different from that of the systems from the seed class, and thereby reveal 
new properties of the deformed systems, which will be discussed in more detail elsewhere. 

1. Main results 

We start with the following simple results that slightly generalize Proposition 2 from [22] and the 
results of [T3|[I2]. The proof is by straightforward computation. 

Proposition 1. Let (M,P) be a Poisson manifold with the Poisson bracket {f,g} = (df,Pdg). 
Consider k functions on M of the form 

H i = HV+ a H$ 1) , i = l,...,k, (1) 

where a is a parameter, and h!- ' and are smooth functions on M. Assume that Hi are func- 
tionally independent for all values of a. 

Suppose that there exists an s G {1, . . . , k} such that Hs 7^ and 

{H s ,H 3 } = (2) 

for all j = 1, . . . , k and for all values of a, and let 

Hi = #f + aH\ 1] = Hl 0) - HV(H® - a)/# s (1) , z = 1, . . . , s - 1, s + 1, . . . , 

H a = Hi 0) +&HP = -{Hi 0) -&)/H^\ (3) 

where a is another parameter. 
Then we have 

{H s ,H J } = (4) 

for all j = 1, . . . , k and for all values of a. 

Corollary 1. Under the assumptions of Proposition 1 suppose that we have {Hi,Hj} = for some 
(fixed) i and j and for all values of a. Then 

{H u H j } = 
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for all values of a. 

The transformation from Hi to Hi is known as a coupling- constant metamorphosis ^3] or as a 
(generalized) Stackel transform P ITo]. 

^From Proposition and Corollary ^ it is immediate that the transformation Q preserves (su- 
per) integrability: if the dynamical system associated with H s is Liouville integrable (so k > n = 
|rankP and H s belongs to a family of n commuting Hamiltonians Hi such that PdHi ^ for all 
i) or superintegrable (i.e., H s is Liouville integrable and k > n), then so is the dynamical system 
associated with H s . 

Note that the relations Q can be inverted: 



Hf ] = Hi- H^H S , i = 1, . . . , s - 1, s + 1, . . . , k, 

h® = -hPh. + &. 



(5) 



Recall that the equations of motion associated with a Hamiltonian H and a Poisson structure P 
on M read (see e.g. j3]) 

dx b /dt H = (X H )\ 6 = 1,..., dimM, (6) 

where x b are local coordinates on M, Xh = PdH is the Hamiltonian vector field associated with H, 
and tn is the corresponding evolution parameter. 

Consider the equations of motion © for H = H s and tu = t and for H = H s and tn = t: 

dx b /dt = (X Hs ) b , 6 = 1, . . . , dimM, (7) 

dx b /di= (X ks ) b , 6 = 1,..., dimM. (8) 

According to [14J we have a reciprocal transformation (see e.g. (201 1221 1231 for more details on such 
transformations) relating the times t and t: 

dt = -H^dt. (9) 

When does © turn © into ©? From © we find that 

d/cft = -H^d/dt, 

and taking into account © and © we see that our question boils down to the following: when does 
the equality 

X Hs = -H^X^ (10) 

hold? 

We have X Hs = PdH s = PdHf ] + aPdH [ P and 

X 6 = PdH s § \^PdH^ + H ^ ~ a 9 PdH?\ 



Plugging this into (|10jl and multiplying the resulting equation by H^, which is nonzero by assump- 
tion, we obtain the following equation: 

(#i 0) + aH^ - &)PdH { p = 0. (11) 

Clearly, PJ holds if and only if either PdH^ = or (#i 0) - 5) = 0. The first possibility 

immediately yields the following result: 

Proposition 2. Under the assumptions of Proposition^ suppose that H^ is a Casimir function 
for P, i.e., PdH^ = 0. Then the transformation (0J) sends the equations of motion {7J) for H s into 
the equations of motion M) for H s . 
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The second possibility is slightly more involved and will be of greater interest to us in the sequel: 

Proposition 3. Under the assumptions of Proposition^ the transformation (GJ) sends the equations 
of motion ffi) for H s restricted onto the level surface H s = a of H s into the equations of motion J3J) 
for H s restricted onto the level surface H s = a of H s . 



Proof. Indeed, if H s = a then + aH^ — a = by On the other hand, the condition H s 
is equivalent to H s = a by (jSJ) • Thus, if H s = a or H s = a then (fTTj) holds, and the result follows. □ 

Note that as the parameters a and a are arbitrary, will transform the equations of motion 
restricted onto any given level surface of H s into the equations of motion (jEJ) restricted onto any 
given level surface of H s . Thus we have a remarkable duality among the deformation parameters and 
the energy (eigen)values that can be readily transferred to the quantum case. 

2. MULTITIME EXTENSION 

Now suppose that all Hi are in involution: 

{H i ,H j } = 0, i,j = l,...,k. (12) 

Then by Corollary so are Hf 

{H t ,H,} = 0, i,j = l,...,k. 
Consider the simultaneous evolutions 

dx b /dt i = {X Hi ) h , b= l,...,dimM, i = l,...,k, (13) 

dx b /di i = (X ili )\ b= l,...,dimM, i = l,...,k, (14) 
and the following extension of (JHJ): 

k 

dt s = - ^H\ 1] dt u ti = U, i = l,2,...,k, (15) 
i=i 

It is straightforward to verify that by virtue of (0), ()12J) . and (|T3*j) we have 

{^°Uf} = 0, {FfUfl^O, (H?\ = (F«) ti , i,j=l,...,k, 

so the transformation (fT3j) for t s is well-defined. This is where we need the commutativity of Hi. 
By virtue of (fT3j) we have 

d rr m d d d tt (d d . . _ . , 

dt s dt s dti dti dt s 

In view of (|T3*|) and (fT4*|) we search for conditions when 

X Hs = -HWx Ss1 X Hi = X 6i -H^X 6s , i = l,2,... ,8-1,8 + 1,. ..,k. (16) 

Plugging (HJ) and (J3J) into (fTSJl yields, after some simplifications, 

(i^i 0) + - a) PdH\ X) = 0, z = 1, . . . , k. 

Hence we have the following generalizations of Propositions 121 and |3J 

Proposition 4. Under the assumptions of Proposition^ suppose that 

{Hi,Hj} = 0, i,j = l,...,k, 

and that H^, i = 1, . . . , k, are Casimir functions for P, i.e., PdH- 1 ^ = 0, i = 1, . . . , k. Then M5\) 
transforms (Q3J) into 
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Proposition 5. Under the assumptions of Proposition^ suppose that 

{H i ,H j } = 0, i,j = l,...,k. 

Then the reciprocal transformation hi 3^ sends the equations of motion (Q3J) restricted onto the level 
surface H s = a into the equations of motion \1J\) restricted onto the level surface H s = a. 

3. Canonical Poisson structure 

Let P be a canonical Poisson structure on M = M. 2n . Then the Hamilton-Jacobi equations for Hi 
and Hi have a common solution, cf. [9]. Namely, we have the following extension of the results of [Oj 
to the Hamiltonians that are not necessarily quadratic in the momenta: 

Proposition 6. Under the assumptions of Proposition^ let M = M. 2n , P be a canonical Poisson 
structure on M, and \, fa, % = 1, . . . , n, be the Darboux coordinates for P, i.e., {Xi, fa} = 5^. Let 
X = (Ai, . . . , A n ) and n = (fa, . . . ,fa n )- 

Let S = S(\,a,E s ,ax,...,a n -i) ! where ai are arbitrary constants, be a complete integral of the 
stationary Hamilton-Jacobi equation for the Hamiltonian H s = H s (a, A, //), 

H s (a, X,dS/dX) = E s . 

If we set E s = at and a = E s then S = S(X,a,E s ,ai,...,a n -i) is a complete integral of the 
stationary Hamilton-Jacobi equation for the Hamiltonian H s = H s (a, X, yb), 

H s {a, A, dS/dX) = E s . 

Moreover, let {Hi, Hj} — 0, i, j — 1, . . . , k, and let 

S = S(X,a,E 1 , . . . , E k , ai, . . . , a„_ fc ) (17) 

where ai are arbitrary constants, be a complete integral for the system of stationary Hamilton-Jacobi 
equations 

Hi(a,X,dS/dX) = Ei, i = l,...,k. 

If we set 

a = E s , E s = a, Ei = E h i = 1, 2, . . . , s - 1, s + 1, . . . , k, 
then S (fT7[ ) is also a complete integral for the system 

Hi(a,X,dS/dX) = E h i = l,...,k. 

As for the equations of motion, we have, in addition to general Propositions El and El a somewhat 
more explicit result: 

Corollary 2. Under the assumptions of Proposition [3| let M = M. 2n , P be a canonical Poisson 
structure on M, and Aj, fa, i — 1, . . . , n be the Darboux coordinates for P, i.e., {Xi,fa} = 5ij. Let 
X = (Ai, . . . , A n ) and /x = (fa, ■■■,fai)- 

Suppose that k = n, dH!> /d/j, = for alii = 1, ... ,n and that Aj can be chosen as local coordinates 
on the Lagrangian submanifold Ne = {(A, /n) G M\Hi(a, X, /x) = Ei, i = 1, . . . , n} (in other words, 
the system Hi(a, X, /j,) = Ei, i = 1, . . . , n can be solved for fi ), and that we have 

E s = a, E s = a, Ei = Ei, i — 1, 2, . . . , s — 1, s + 1, . . . , n. (18) 

Then the reciprocal transformation (Q3J) turns the system 

dX/dU = (dHi/d/j,)\N E , i = l,...,n, (19) 

into 

dX/dU = (dHi/dn)\fj_, z = 1, ... ,n, (20) 
where N E = {(A, fi) e M\Hi(a, A, fx) — Ei, i — 1, . . . ,n}. 



6 MACIEJ BLASZAK AND ARTUR SERGYEYEV 

For instance, if we have 

Hi = ^0, Gi(X)n) + Vi(X) +aWi(X), i = l,...,n, (21) 

where (•, •) stands for the scalar product in R n and Gi(X) are n x n matrices, then the system ()19|) 
reads 

dX/dti = Gi{X)M, (22) 
where /i. = M(X, a, E\, . . . , E n ) is a general solution of the system H^a, X, fi) = E iy i — 1, . . . , n. 
If we eliminate M from (}2*2*|) then we obtain the dispersionless Killing systems (cf. jH]) 

X ti =G i (G s )- 1 X ts , i = l,2,...,s-l,s + l,...,n, (23) 
and the reciprocal transformation (fTo^) . which in our case reads 



dt s = - } u Wi(X)dti, U = t { , i ^ s, 



i=l 

turns (j2H| into 

X k =G t ,(G s )- 1 X- ts , i = l,2,... ,8-1,8 + 1,. ..,n, (24) 

where G s = — G s /W s and Gi = Gi — WiG s /W s , i = 1, 2, . . . , s — 1, s + 1, . . . , n, are related to the 
Hamiltonians 

Hi = i(/x, Gi(X)fi) + Vi(X) +aWi(X), i = 1, . . . , n. (25) 

4. Solving the reduced equations of motion 

Under the assumptions of Corollary |21 we can apply Proposition |B] in order to obtain the solutions 
of equations of motion (fTU|) and (|2*Uj) as follows: 

Corollary 3. Under the assumptions of Corollary^ suppose that 

S = S(\a,E 1 ,...,E n ) (26) 

is a complete integral for the system of stationary Hamilton- J acobi equations 

H(a, X,dS/dX) = Ei, i = l,...,n. 

Then a general solution of M9\) for i = r can be written in implicit form as 

OS/dEj = 5 jr t r + bj, j = 1, . . . , n, (27) 

where bj are arbitrary constants, and by virtue of < f7^)) a general solution of H?J) for i = r can be 
written in implicit form as 

dS/dEj = 5 jr t r + bj, j = 1, . . . , n. (28) 

Comparing (|2*7j) and (J2~%j) and using (|TS|) we readily see that, in perfect agreement with tj = ti 
for i ^ s, but t s = dS/ dE s — b s = dS/ da — b s while t s = dS/ dE s — b s = dS/ da — b s , so this approach 
does not yield an explicit formula expressing i s as a function of A, fx, and ti. 

In order to find a complete integral (|26|) we can use separation of variables as follows (see e.g. 
j7] and references therein). Under the assumptions of Corollary 01 suppose that A i; /ij, i = 1, . . . ,n 
are separation coordinates for the Hamiltonians Hi, i = 1, . . . , n, that is, the system of equations 
Hi(a, X, n) = Ei, i — 1, . . . , n, is equivalent to the following one: 

(fi(Xi, fj,i, a, E x , . . . , E n ) = 0, i = l,...,n, (29) 

that is, the separation relations on the Lagrangian submanifold Ne- 
Consider the system of stationary Hamilton-Jacobi equations for Hi 

Hi(a,X,dS/dX) = E u i = l,...,n. (30) 
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By the above, (jSOjl is equivalent to the system 

ip i (\ i ,dS/d\ i ,a,E 1 , . . . ,E n ) = 0, i = l,...,n. (31) 

Suppose that (f2~H|) can be solved for fa, i — 1, . . . , n: 

Hi = Mi(Xi, a, E h . . . , E n ), i = 1, . . . , n. 

Then there exists a complete integral of (|3*Tj) . and hence of of the form (cf. e.g. [2] and 

references therein) 

n „ 

1=1 ^ 

and general solutions for (fTH|) and (f2T)|) can be found using the method of Corollary El 
The formulas (f2~Tj) take the form 

n „ 

/ (dMi(\i, a,E l7 ..., E n )/dEj)d\i = 5 jr t r + b j} j = 1, . . . , n. 
i=i J 

For r = s we have 

n „ 

£ s + 6 s = dS/dE s = dS/da = J2 ( dM i( X ^ a > E i> • • • > E n )/da)d\i. 

8=1 ^ 

5. Deformations of seed systems 

Under the assumptions of Corollary 2, suppose that Aj, /ij, i = 1, . . . , n are separation coordinates 
for the Hamiltonians Hi, i — 1, . . . ,n, then the Lagrangian submanifold iVs is defined by n separation 
relations (jHU). Further assume that all functions (pi are identical, 

ifi = (fi(\i,Hi, a,Ex,... ,E n ), i = l,...,n. (33) 

Then the relations (J2U)) mean that the points \i, fa belong to the separation curve 

<p(\,fa a,Ei,..., E n ) = (34) 

for alH = 1, . . . , n. 
If the relations 

(p(\i, fa, a, Hi,..., H n ) = 0, i = 1, . . . , n, 

uniquely determine the Hamiltonians Hi for i = 1, . . . , n, then for the sake of brevity we shall say 
that Hi for i = 1, . . . , n have the separation curve 

ip(\,H,a,Hi, . . .,H n ) = 0. 

Fixing values of all Hamiltonians Hi = E iy i = 1, . . . ,n, picks a particular Lagrangian submanifold 
from the Lagrangian foliation. Setting a = in the above formulas we see that Aj, fa, i = 1, . . . , n 
are separation coordinates for the Hamiltonians Hf \ i — 1, . . . , n, as well, so the system of equations 
if f (A, ^t) = EJj, z = 1, . . . , n, is equivalent to 

ip o (\i,fa,Ei,...,E n ) = 0, i = l,...,n, (35) 

where (p Q (\i, fa, E x , . . . , E n ) = (p(\ Ei, . . . , E n )\ a=0 . 

In what follows we shall restrict ourselves to the separable systems whose separation curves (po = 
read 

tfJV 1 + #< 0) A* + • • • + Iff A^" = V(A, //) (36) 

where n + A; — 1 = (3i > fa > ■■■ > Pn = 0, k 6 N and -?/>(A,/i) is a smooth function. Each class of 
systems (J31j|) is labeled by a sequence (Pi, (3 n ) while a particular system from a class is given by a 
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particular choice of //). In particular, the choice ^(A, /i) = |/(A)/i 2 + 7(A) yields the well-known 
classical Stackel systems. 

For k = there is only one class given by 

H^X™- 1 + Hf ] \ n - 2 + ...+H® = ip(X, y) (37) 

which is precisely the Benenti class of Stackel systems jT| |2] if ip(X, y) = |/(A)/i 2 + 7(A). All these 
systems separate in the same set of coordinates (Aj,//j) by construction. We shall refer below to the 
systems with the separation curve ()37|) as to the systems from the seed class. 

In [7j it was shown that an arbitrary class of the systems with the separation curve (|3l)j) is obtained 
via the so-called /c-hole deformation from the seed class (|37p. 

Below we demonstrate that an arbitrary /c-hole deformation is nothing but a sequence of k Stackel 
transforms (J3J), and hence all separable classes (J31j|) are Stackel-equivalent to the seed class in the 
sense of j2] . In order to do this we first introduce an alternative notation for different classes (J36j) . 
which is more convenient for further considerations as well as for the bi-Hamiltonian extension. 

s 

We shall call a polynomial of the form ^ a m+ jX n+k ~^ m+ ^ a sub-chain of length s if a& 7^ for 

3=1 

k — m + 1, m + s, and a string of holes of length s if = for k — m + 1, m + s. 

n+k 

Now consider a polynomial Yl CLi.X n+h ~ l , where precisely k of the coefficients cij equal zero, but 

i=i 

cii 7^ and a n+ k 7^ 0. Denote n non- vanishing coefficients by {h[°\ Hn^), where 7^ a\ = 
and 7^ a n+ k = Hn^ by assumption. In what follows we also assume that if 7^ at = H$ and 
^ ai = HP and k > I then k' > V . 

It is immediate that any class of separation curves (J31j|) is uniquely determined by a sequence of sub- 

1 1-1 

chains and strings of holes (72,1, m 1; n 2 , m 2 , mi_i,ni), m, = n,Y2 m i — k, where is the length 

i=l i=l 

of i-th sub-chain and m,j is the length of i-th string of holes. The separation curve corresponding to 

1 1-1 

a sequence (ni,mi, n 2 , m 2 , m;_i, ni), Hi = n,J2 rrii = k, reads 

i=l i=l 

± HfX^-> + £ "f ff<°> , A ««-£<«.«-.+"-.-) + £ A „ 1 - j = v , (Aj „ } . (3g) 

For fc = we have only one chain (the seed class). For k = 1 there are (n — 1) different classes 
consisting of two sub-chains (72,1,1,712), % + n 2 = n, separated by one hole. For k = 2 we have 
\{n — l){n — 2) different classes, where (n — 1) of these classes consist of two sub-chains separated 
by a two-hole string (72,1,2,712), rii + n 2 = n, while the remaining cases consist of three sub-chains 
(72,1, 1, n 2 , 1, n 3 ), ni + n 2 + n 3 = n, separated by single holes, and so on. 

Now define the transformation from the /c-hole case to the (k + l)-hole one. Without loss of 
generality we can restrict ourselves to considering the following subcases only: 

(i) (m, mi, n 2 , m 2 , m/_i, n[) -> (72,1, mi, n 2 , m 2 , m z _i + 1, m), m = n\, 

(ii) (ni, mi, n 2 , m 2 , mj_i, nj) -> (ni, mi, n 2 , m 2 , m z _i, n z , 1, n z+1 ), m + m+i = rij. 

In the case (i) the number of sub-chains and their lengths are preserved, while the length of the 
last string of holes is increased by one. In the case (ii) the last sub-chain is split into two sub-chains 
by inserting an additional hole. Notice that we can reach an arbitrary s-hole deformation in a unique 
way from the seed class (n) by applying the above recursion step s times. 

Passing from the fc-hole deformation to the (k + l)-hole one means, according to our recursion, 
that for the separation curve we have 

rr(0) ^n+k-i fr(0) ■> n +k-i+l ■ _ i 

H i x —> ti i x , % — 1, . . . ,m+i 
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fffV^tffVs i> n , + i. 
If rii = and rii + i = n\ then we have the case (i) while for n; +1 < n\ we have the case (ii). 

For the sake of convenience we now formally merge the cases (i) and (ii) into a single transformation 

(ni, m 1 ,n 2 , m 2 , ■ ■ ■ , raz-i, nj) -> (ni, mi, n 2 , m 2 , . • . , raj-i, nj, 1, n J+ i), m + n i+1 = n[, 

where n\ = for the case (i) and 1 < n\ < n[ — 1 for the case (ii). Here n\ = corresponds to a void 
sub-chain (sub-chain of zero length). 

Proposition 7. Consider two n-tuples of Hamiltonians {i^ } and {f^ ' 1 } iae separation curves 
of the form 

f tff A^- + g V # (0) r A n+fc - -1^-.^,-. ) + £ ^ m _, = ^ ^ 

(39) 

corresponding to the k-hole deformation of the seed class and 



-(ri\ , , ■ ni - fri\ , ■ l-mi+l-r _ n+k+1- T, (n; +1 _ s +TO;_„) 

£ ffj°)A»+*+i-J + £ rfi+i A ~ n,+1 ~ ; ' + E E # ( } r A <=1 

i=l j=l r=2 i=l /^q\ 

+ E^i 1+J A ni - J = ^(A,/i), 

corresponding to the (k + l)-hole deformation of the seed class Jff7| ), respectively. Here n\ = and 
n^ + i = /or tae case fa) and 1 < n; < n^ — 1 and nj = n; + i + n; for the case (ii); mi = 1 in both cases. 

Then the Stackel transform from {H^} to {H^} reads 



rr(0) = 1 H-(O) 

1 ..fa^itl nu 



T/ (n+fc) 

y(n+fe) 



flf = - f5^»£>, < = 2- (41) 



T/ (n+fe) "i+i 
rr(n+k) 

fr _ o-(O) u-(o) • 



where vj; n+k ^ are separable potentials defined by the relation [Z| 

(42) 

that must hold for A = Aj, i — 1, . . . , n. 

Proof. In order to compare the Hamiltonians we should reduce the separation curve Pjl for Hf ] to 
that for -f/j-°' ) . To this end we get rid of the highest monomial H^\ n+k in f40|) by expressing \ n+k 
from (|4"2*|) . Then comparing coefficients of the separation curves for and yields 

H^ = H^-Vr k) Sr\ i = l,...,n l+1 -l, 

Hit = -y^ + \ k) Hf\ 

Hf = H®-Vp**>H®, i = n l+1 + l,...,n 
and (jUJ) follows. □ 
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Notice that after the remuneration of Hamiltonians Hf , 



h[ 0) = H® , 
1 n i+i > 

Hl 0) = Ht\, i = 2,...,n l+1 , (43) 
Hf^Hf\ i>n l+1 , 

we deal with a particular case of Proposition [H where = V^ n+k \ a = 0, s = ni + \. Thus, the 
reciprocal transformation 

n 

di n l+1 =~Y1 Vt^tij' *< = M ^+1' ( 44 ) 



transforms the equations of motion for Hamiltonians Hi — H- + ai^ of the /c-hole deformation 



restricted onto the level surface H ni+1 = into the equations of motion for Hamiltonians H^ of the 
(k + l)-hole deformation restricted onto the level surface Hni +1 = h[ ^ = a. 
We can readily extend to Hi using (J3J) and the results of 0. Namely, let 

fH 1 ) _ _ _ y(n+k-m +1 ) 

Then the separation curves for ifj and if; are of the form 

a\ n+k + iJ ni+1 A n+fc ~ n ' +1 + ^ H t X^ = i/;(X,fi), 

H n , l+1 X n+k + aX n+k ~ ni+1 + =^{X,fi), 

and the reciprocal transformation (}4*4*|) transforms the equations of motion for Hamiltonians Hi 
restricted onto the level surface H ni+1 = a into the equations of motion for Hamiltonians Hi restricted 
onto the level surface H ni+1 = a. 

As the whole procedure is recursive, it means that the n-tuple of integrable Hamiltonian dynamical 
systems described by the separation curve (J36|) . restricted onto a given Lagrangian submanifold, are 
related to an n-tuple of Hamiltonian dynamical systems from the seed class (also restricted onto an 
appropriate Lagrangian submanifold) via the sequence of reciprocal transformations. 



6. Example 

As a simple illustration of the above results, consider the Henon-Heiles system on a four- dimensional 
phase space with the coordinates (pi,P2> <Zi> (Za) and canonical symplectic structure. The correspond- 
ing Hamiltonian 

Hi = \v\ + \v\ + Qi + ^QiQl ~ Oiqi, 



is in involution with 



u 1 2 ^12, ^ 4,-^22 1 2 
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Consider the following one-hole deformation (Stackel transform) of Hi and H 2 : 

Hi = -—Pi + -— P 2 + <?i + 7"<? 2 - «— 



qi zqi i q l 

?2 2 1 <?2 2 1 4 i ~ 1 ?2 

gi Sg 1 16 4g 2 

The corresponding reciprocal transformation 

1 



^2 = 2^2 - 2^2 - g-Pi - g-^2 - 7^2 + ^-T- 



dti = q x dti + ^qjdt 2 , t 2 = t 2 

defines the map between equations of motion restricted to the respective Lagrangian submanifolds, 
N E = {Hi = a,H 2 = E} and = {Hi = a,H 2 = E}. 

The separation coordinates (A, /i) are related to p's and g's by the formulas 



qi = Xi + A 2 , g 2 = 2V-A1A2, 

Ai/ii A 2 /i2 / — > — < / A*i A*2 

Pi = 1 — + \ i - : = V -AiA 2 T r- + 



Ai — A 2 A2 — Ai \Ai — A2 A2 — Ai 

are common for Hi, H 2 and H 2 and the respective separation curves read 

a A 2 + Ai2i + H 2 = \x^ 2 + A 4 , 

H1X 2 + aX + H 2 = ^A/i 2 + A 4 . 
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